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A transition between energy levels at an avoided crossing is known as a Landau-Zcner transition. 
When a two-level system (TLS) is subject to periodic driving with sufficiently large amplitude, 
a sequence of transitions occurs. The phase accumulated between transitions (commonly known 
as the Stiickelberg phase) may result in constructive or destructive interference. Accordingly, the 
physical observables of the system exhibit periodic dependence on the various system parame- 
ters. This phenomenon is often referred to as Landau-Zener-Stiickelberg (LZS) interferometry. 
Phenomena related to LZS interferometry occur in a variety of physical systems. In particular, 
recent experiments on LZS interferometry in superconducting TLSs (qubits) have demonstrated 
the potential for using this kind of interferometry as an effective tool for obtaining the parameters 
Q 1 characterizing the TLS as well as its interaction with the control fields and with the environment. 

q ' Furthermore, strong driving could allow for fast and reliable control of the quantum system. Here 

I ' we review recent experimental results on LZS interferometry, and we present related theory. 

5-H ' 

C/3 



Contents 



I. Introduction 2 

' O ' A. Strongly driven TLSs: a brief historical review 2 

1. Atomic physics 2 

O , 2. Other systems 3 

O i 3. Theory of LZS interferometry 2 

B. Superconducting qubits 2 

C. Fano and Fabry-Perot interferometry using superconducting qubits 4 

D. Hamiltonian and bases 4 

> , 

I"""""'" , II. Theory: adiabatic-impulse model 6 

A. Adiabatic evolution 7 
■ B. Single passage: Landau-Zener transition 7 

C. Double passage: Stiickelberg phase £ 

D. Multiple passage £ 

1. Slow-passage limit 11 

2. Fast-passage limit 12 

E. Decoherence [13 
' F. "Breakdown" of the adiabatic theorem 17 

G. Reversal of stimulated emission llg 

k>( i III. Recent experiments with strongly driven superconducting quantum circuits 1£ 

I^N , A. Multiphoton transitions in superconducting qubits 1£ 

B. LZS interferometry in superconducting qubits 1£ 

C. Multiphoton resonances in multi-level systems [21 

IV. Conclusions [13 

Acknowledgments [2l] 

A. Solution of the Landau-Zener problem 21 

1. Adiabatic wave function 21 

2. Non-adiabatic transition (evolution matrix for the Landau-Zener transition) [22 



'Electronic address: sshcvchcnkoOilt. kharkov.ua 



2 



B. Evolution of a periodically driven two-level system 



C. Rotating-wave approximation 

1. Hamiltonian in the RWA-form 

2. Solving the Schrodinger equation in the absence of relaxation 

3. Solving the Bloch equations with relaxation 

D. Floquet theory 

E. Dressed-state picture: quantized driving field 
References 



2fr 

26 

21 

27 
301 

:3d 



I. INTRODUCTION 

A. Strongly driven TLSs: a brief historical review 

The model of a quantum two- level system (TLS), with only two relevant quantum states, is used to describe a 
variety of physical systems. At first, it was used in relation to spins and atomic collisions, which describe natural 
microscopic systems. It has been used to describe various other naturally occurring settings ever since. Furthermore, 
some artificial mesoscopic systems that were realized recently in semiconductor quantum dots and superconducting 
circuits can operate as effective TLSs. These solid-state TLSs have received increasing attention recently for two 
main reasons: they exhibit fundamental quantum phenomena on macroscopic scales, and they are considered possible 
candidates to operate as quantum bits (qubits) in future quantum information processing devices. 

It is quite common that the two energy levels of a quantum TLS exhibit an avoided level crossing or anticrossing 
as some external parameter is varied, as shown in Fig. [T] The physical properties of the two energy eigenstates 
are typically exchanged when going from one side of the avoided crossing to the other side. If the external control 
parameter is varied in time such that the system traverses the avoided crossing region, a non-adiabatic transition 
between the two energy levels can occur. The transitio n probability is usuall y named aft er Landau and Zener (LZ) 
in recognition of their pioneering work on this subject, ( Landaul Il932tj) and ( Zenerl . [l932f) . 

If a LZ transition is one step in a sequence of coherent processes occurring during the dynamics, not only the change 
in the occupation probabilities but also the change in the relative phase between the quantum states is relevant. In 
particular, in this review we shall focus on the situation where the control parameter is varied periodically and 
strongly, such that the system keeps going back and forth across the avoided crossing reg ion. In this case , the physical 
obser vables exhibit periodic dependence on the phase acquired between the transitions ( Landaul Il932al Istuckelbergl . 
Il932h . This periodicity, which is the basis of Landau-Zener-Stuckelberg (LZS) interferometry, provides a useful tool 
that allows the characterization of the parameters defining the quantum TLS and its interaction with the control fields 
and the environment. We not e here that the formula for the single-passage transition probability can also be found 



in the paper (IMaiorana, 1932), and it can alte rnatively be named the Landau-Zener-Stiickelberg-Majorana formula, 
as proposed in ( Di Giacomo and Nikitinl . |2005|) . 



1. Atomic physics 

The problem of strongly driven TLSs appears in a number of contexts in atomic physics. In atomic-collision 
problems, the electronic state of the atom pair can experience a LZ transition from the ground states of the incoming 
atoms into a hybridized intermediate state and t hen back to the electronic ground states of the outgoing atoms. In this 
case the avoided crossing region is passed twice ( Child! . Il974t iNikitin and Umanskiil . 11984 ) , and the Stuckelberg phase 
typic ally has large variations and hence its interferometric consequences are washed out ( Landaul Il932al IStuckelberd . 
1 19321) . 

An atom in an intense laser field can also experience LZ transitions, but in contrast to atomic-collis ion problems, 
the atom is driven periodically in time and multiple LZ transitions occur (jDelone and Krainovl 119851 ). In the case 
of an atom in a laser field, (regardless of the field intensity) if the atom is driven at a frequency that matches 
its energy-level splitting, it can be resonantly excited to an upper le vel. Then, at low driving amplitudes, the 
level occupation probabilities experience Rabi oscillations (jRabil . fi"937| ). When the driving amplitude is increased, 
multiphoton processes beco me relevant. T hese processes occur when the energy-level splitting matches the energy of an 
integer number of phot ons (IShirlevl 19651). These m ultiphoton pr oc esses have tra di tionally been desc ri bed by making 
use o f Floquet theory [dAutler and Townesl Il955h. dRitus). Il967l). dSambd. I l973l) (|Zerdovichl . Il973l ). (jBarone et all 
I1977D . (|Aravind and Hirschfelderl . Il984l) ; for a review see (|Chu and Tel'novi . 12004) ] . In addition to controlling atomic 
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states, resonant driving of atom s has also been studied as a means for findi ng atomic parameters, and the parameters 
of quantum systems in general ( Coffey et all Il969t iHenrv and Land . 11977ft . 



2. Other systems 

In addition to atomic systems, LZS interference has been studied in a number of other systems. For example, 
periodic dependence and time-domain osc illations related to LZS in terference have been predicted in the electronic 
transport in semiconductor superlattices ( Rotvig et all Il995l Il996() . Stiickelberg oscillations were also studied in 
super conducting quantum point contacts biased with a dc voltage and exposed to periodic driving ( Gorelik et all 
1998). This effect was proposed as a sensitive tool for microwave-photon d etection. Another examp le of interference in 



microscopic systems aris es in the study of interacting magnetic molec ules [ (jWernsdorfer et al. , 2000), (|GaranirJ . l200l ]. 
molecular nanomag nets ( Calero et all . 120051 : iFoldi et all . 120071 l2008f ) , and single nitrogen vacancy centers in diamond 



(Fuc hs et flE 1200a) . Interference effects have also been pred icted to take place in the inter la yer tunneling in quasi-one- 
dime nsional layered materials, such as organic conductors [( Cooper and Yakovenkol . l2006ft . dBaneriee and Yakovenkol 
12008ft ]. and in the tunneling of a periodically driven two-mode Bose-Einstein condensate (Zh ang et al\. | 2008allbft . 
Recen tly, LZ tunneling and Stiickelberg oscillations were measured in a gas of ultracold molec ules (jMark et ali 
l2007al lbf) and in the dipole-dipole interaction between Rydberg atoms (|van Ditzhuiizen et all . 1200a). 



Avo ided level cr o ssings can also be used in the study of second-order phase transition problems ( Damski and Zurekl . 



200f> ; iDziarmagai l2009ft. Periodic driving around a phase transition can then result in interference fringes 
( Mu kheriee and D utta , 2009). The LZS model has also been use d to describe the singlet-tri pl et transition for th e 
state of two ele ct rons in a sem i condu ctor double quantum dot [( Ribeiro and Burkardl . [2009ft . ( Petta et all . l20ldft . 



(Burkar dl. l2010ft . (Rib eiro et al ., 2010)1 as well as the tr ansport of photons between two halves of an optical cavity 
divided by a dielectric membrane ( Heinrich et all l201dft . Recent technological advances have also allowed the real- 
ization of the strong-driving regi m e in superconducting q u bits, where such inte rf erence has been observ e d in severa l 
experi ments fdOliver et aU l2005ft. (|Sillanpaa et aLl . l2006ft . (|Wilson et aLl . l2007ft . ([izmalkov et aUl2008ft . (|Sun et all . 
2009), (jLaHave et all 120091 )]. Such superconducting qubit systems will be discussed in more detail below. 



3. Theory of LZS interferometry 



The theory of driven qu an tum TLSs has been st u died in numerous pap e rs [for reviews se e (iLeggett et afl .ll987ft. 
(|Grifoni and Hanggil [l998ft . (|Thorwart et aLl . l200l[ ). (|Vitanov et a /ll200l ft. dZhu et all. [2001ft . jFrascal l2003ftl . Some 



1991ft . (|Thorwart et all 



of th e recent work has fo cu sed on the impact o f the environment, e.g. [flAo and Rammed . 

l2000ft . (|Wubs et aLl . l2006ft . (|Zueco et ^.1 , 120081) . (|Nalbach and Thorwaru. 12009ft ] . In particular, some papers pointed 
out that if the relaxation in the s ystem is small enough, interf er ence between multiple LZ transitions becomes r elevant 
for the system's dynamics, e.g. ( Shimshoni and Gefenl . Il99"lft . ( Saito and Kavanumal 120021 : ISaito et al. 1. 12007ft : when 
the avoided crossing is passed several time s, interference can take place and the dynamics is sensitive to the phase 
acquired between successive crossings, e.g. (|Kavanumalll994ft . 

LZS interference has also been recently proposed as a possible approach to achieve pe rfect population transfers 
and to implement quant um gates for qua nt um-control and q uantum- computing purposes [( Teranishi and Nakamural 
11998ft . (|Zhu et all l200lft . (|Gaitanl . 12003ft . (|Zwanziger et all . l2003ft . (|Nagava et a ll. 12007ft. dLi et all l2009ftl [related 
work o n combining adiab a tic LZ crossings wi th coherent dyna mics can be found in (|Goswamil . 12003ft . (IVitanov et all 
12003ft . dKral et ad l2007h . (lOh et all . [2008ft . dWei et all [2008ft ], 



The purpose of the present review is to present the theory of strongly driven TLSs and LZS interferometry and 
recent experimental results on the subject, particularly in relation to superconducting qubits. 



B. Superconducting qubits 

Superconducting qubits are Josephson-j unction-based circuits that can behave as effective two-level systems. There 
are several types of superconducting qubits, differing mostly by t heir topology and by t h eir physical parameters. The y 
are generally called p ha se, charge, or flux qubits . For re views see [dMak hlin et al. [ 1 20011) , toevoret and Martinisl.l2004ft . 
(|You and Noril , l2005ft . (jWendin and Shumeikol 12007ft . (|Zagoskin and Blaisl . 120071 ), and (jdarke and Wilhelml , l2008ft ] . 

Superconducting qubits have several properties that make them unique for studying quantum effects. The qubits 
have mesoscopic size, which means that they allow the observation of quantum effects on a macroscopic scale. They 
can be integrated in electrical circuits, which provides different ways of controlling and probing their states. Their 
parameters are adjustable: they can be tuned by changing the applied bias current, gate voltage, or magnetic flux. 
The qubit's state can be probed by measuring the charge or current induced in it. These unique properties of 
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superconducting qubits allow the realization of the strong-driving regime, and they have allowed the realization of 
several experiments on LZS interferomctry in the past few years. 



C. Fano and Fabry-Perot interferometry using superconducting qubits 

Superconducting qubits can be used for LZ S interferometry, b ut also to study othe r types of interferometry, in- 
cluding Fano and Fabry-Perot interferometry [( Zhou et all l2008al lbl). ( Liao et all . l2010t )]. These latter devices can be 



implemented using quasi-one-dimensional open systems where photons are injected from the left and move towards 
the right side of the device. Along the way, the photons interact with either one or two qubits acting as tunable 
mirrors, controlled by changing the applied electric and/or magnetic fields on the qubits. These qubits, working as 
tunable mirrors, can change the reflection and transmission coefficients of the photons confined in waveg uides. 



Let us first consider the case of a single superconducting qubit interacting with an incoming photon ([Zhou et all 
2008a). When the energy of the incoming photon matches the energy spacing of the qubit, the photon is reflected, 
otherwise it is transmitted. This type of single-photon switch exhibits Breit-Wigner scattering: now in one dimension 
instead of the standard three dimensional case for natural atoms. This Breit-Wigner scattering produces a symmetric 
Lorentzian peak in the reflection coefficient, versus frequency, of the photon. This situation occurs when the dispersion 
relation of the incoming photon is linear, as in a transmission line resonator, acting as a "rail" guiding the motion 
of the photons. When the photon dispersion relation is nonlinear, for long-wavelength photons propagating in a 
quasi-one-dimensional array of coupled cavities, the reflection coefficient exhibits an asymmetric Fano line shape, due 
to the interference between the continuous mode of the incoming photon and the discrete energy levels of the qubit. 
Thus, for single-photon transport in a one-dimensional waveguide, the photons can be partially or totally reflected by 
a controllable two-level system which can act as a tunable mirror. 

It is known that the Fabry-Perot interferometer, which consists of two highly reflecting planar mirrors, provides the 
simplest cavity. It is then natural to ask the question: "is it possible to construct a quantum version of a Fabry-Perot 
interferometer?" Namely, to bui ld a resonator, in a one-dimensional waveguide, with two tunable-mirrors made of 



quantum scatterers. Reference (jZhou et all l2008bf ) focused on this question and studied quantum analogs of the 



Fabry-Perot interferometer. They found that two separate two-level systems intera cting with photon s in a waveguide 
can also create, between them, single-photon quasi-bound states. A recent work ( Liao et all . |2010| ) studied a more 



complex form of interferometry, with the goal of controlling the transport of single photons by tuning the frequency 
of either one or two cavities in an array of coupled cavities guiding the ID motion of photons. In some regimes, the 
photons can be lo calized around the sc atterers, which ac t as impurities, producing isolated states in the gap of the 
energy spectrum [ (jZhou et aUl2008bl) . (|Liao et aZ.l . l201Cj) ]. 



D. Hamiltonian and bases 

The main subject of this review is a strongly and periodically driven TLS. A quantum TLS with energy bias e and 
tunnelling amplitude A is described by the Hamiltonian 



H(t) = ~a x - £ Ma z (1) 

in terms of the Pauli matrices a x z (H = 1 is assumed throughout). One generally thinks of A as being fixed by the 
system properties, and e as being a tunable control parameter; hence £ is a function of time t and A is not. We 
assume the monochromatic time-dependent bias 

e(t) = e Q + Asinut, (2) 

with amplitude A, frequency oj and offset eq. 

Note that this problem with the time-dependent component along the z-axis can be related to the one with the 
time-dependent component along the x-axis by a ir/2- rotation around the y-axis: 

H'(t) = e-^H(t)e^°y = -^a z + ^-a x . (3) 

The latter Hamiltonian is typical for problems such as an atom in a laser field. It should also be noted that having 
different signs in the Hamiltonian than the ones given above does not cause any nontrivial change in the results. 
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The instantaneous eigenvalues of H(t) depend on the bias e(t) as follows: 



E ± (t) = ±-n(t), 



The e-dependence of the energy levels is shown in Fig. [TJ In particular, it shows an avoided level crossing at e = 0. 



(4) 
(5) 




FIG. 1 (Color online) (a) Energy levels E versus the bias e. The two solid curves (red and blue) represent the adiabatic energy 
levels, E±, which display avoided crossing with energy splitting A. The dashed lines show the crossing diabatic energy levels 
Ef i, corresponding to the diabatic states tpf and <p±. (b) The bias e represents the driving signal, and it oscillates between 
£min = £o — A and £max = £o + A with a sinusoidal time dependence: e(t) = £o + Asinwt. 



Although the TLS is a rather simple-looking model, concrete calculations on it can present significant difficulties. 
Specifically, the Schrddinger equation can be written as a second- order differential equation w ith periodic coefficients, 
the Hill equation, which is not solvable in analytic closed form (jGrifoni and Hanggil .[l998). Nevertheless, different 
theoretical approaches can be used to obtain approximate analytical results in various parameter regions. 

As we shall explain in detail below, one key relation for purposes of deciding the suitability of a given theoretical 
approach is that be tween the driving freq uency to and the minimal energ y splitting A. In superconducting qubits, 
both limits, u < A (jSillanpaa et aZll2006h and w » A (jOliver et all l2005l i . have been realized in recent experiments. 
Another key relation concerns the velocity of passing the avoided level region, which is of the order of Aoj (we also 
introduce the adiabaticity parameter 5 later in Eqs. (|19H20j) ). We will distinguish between the slow driving regime, 
Aoj < A 2 , and the fast driving regime, Aoj 3> A 2 . The natural representations for the two regimes are formulated 
in the adiabatic and diabatic bases, respectively. The diabatic basis {ytiVj.} is formed with the eigenstates of a z : 
cr z tp-f = </5f , u z ip^ = — ipi (note that these states would be the eigenstates of the Hamiltonian if A vanished). The 
adiabatic basis consists of the instantaneous eigenstates of the time-dependent Hamiltonian: H(t)tp± (t) — E± (t)tp± (t) . 
The two bases are related to each other by the time-dependent coefficients /3±(t) (see Appendix A for details): 



<P±(t) = p T (t) n T P±(t)n, 

(3±(t) = 



'n(i)±e(t) 



2Q(t) 



(6) 
(7) 



Any given wave function of the quantum TLS can be decomposed in either one of these two bases; the coefficients in 
the quantum superposition are the probability amplitudes for the respective states. 
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II. THEORY: ADIABATIC-IM PULSE MODEL 



There are a number of theoretical approaches that can be used to study the problem of a strongly driven TLS. 
Among these approaches, the adiabatic-impulse model, which we shall explain in detail below (and in Appendices A 
and B), is perhaps the most intuitive whenever the TLS is repeatedly driven through the avoided crossing. We shall 
therefore put special emphasis on this theoretical approach in this review. Three alternative theoretical approaches, 
which in certain parameter regimes can be well suited for the study of this problem, will be discussed in Appendices 
C (rotating- wave approximation), D (Floquet theory) and E (dressed-state picture). 

The adiabatic energy levels E± introduced in Sec. I.C have a minimum distance of A, and this minimum distance 
is realized at times t\$ + 2irn/uj, where uti = arcsin(— Eq/A) and ujt2 = n — ojt\, see Fig. [5] here n is an integer. 
Because the energy eigenstates, i.e. the states of the adiabatic basis, change rather rapidly around the avoided-crossing 
region and are approximately constant far to the right or far to the left of that region, one would intuitively expect 
that the system will evolve almost adiabatically far from the points of avoided crossing and the evolution becomes 
nonadiabatic in the vicinity of these points. One could therefore make the approximation that beyond a certain 
boundary the evolution is completely adiabatic. However, it is convenient not to introduce an arbitrary boundary 
between the adiabatic and nonadiabatic regions and instead consider the evolution to be adiabatic everywhere except 
at the points of minimum energy separation, where the system experiences a sudden mixing in the populations 
of the two energy levels. Naturally, these sudden transitions do not occur in reality. However, provided that the 
system follows a linear ramp traversing the avoided crossing region, one can still obtain an accurate description of 
the net result of the adiabatic and nonadiabatic parts of the evolution. This discretized picture only simplifies the 
algebra that arises in the relevant calculations. The calculation presented in Appendix A can be used to obtain the 
appropriate assignment for the adiabatic phases and the mixing matrices at the crossing p oints in o rder to obtain the 
correct express ions describing the system's dynamics. Following ( Damski and Zurekl 120061 ) (see also ( Zurekl . [i996|) and 
( Damskil . 120051 )) we call this picture the adiabatic-impulse approximation to emphasize the two-stage character of the 
model. In particular, this name emphasizes that the non-adiabatic transitions are described as instantaneous, which 
is just a convenient d escription of the continuou s dynamics; see also original artic le s, where the adiabatic-impuls e 
theory was developed (iDelone and Krainov . 19851). dAverbukh and Perel'manl . fl985l) . ( Vitanov and Garrawavl [l996l) . 
( Garrawav and Vitanov , 1997 ). ( Damski and Zurekl 20061) ]. 
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FIG. 2 (Color online) Time evolution of the energy levels during one period. The time-dependent adiabatic energy levels 
define a two-stage evolution: transitions at the non-adiabatic regions, described by the evolution matrix N, and the adiabatic 
evolution, described by the matrices Ui t 2 = exp(— iCi,20" z ). The acquired phases £i,2 have a geometrical interpretation: they 
are equal to the area under the curves, shown by the yellow and blue regions. The diabatic energy levels, ±e(t)/2, are shown 
by the dashed lines. 
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A. Adiabatic evolution 

The wave function ip(t) that describes the quantum state of the TLS as a function of time can be decomposed in 
the adiabatic basis (see Appendix A): 

m = 5>±(%±(t) = ^c ± (t)<Mt)exp{^(c + J)}, ( 8 ) 
± ± 

C = \J (l(t)dt. (9) 

Since in the adiabatic approximation c± are time-independent coefficients between crossings, the adiabatic evolution 
from t = ti to t = ti (assuming that no LZ transitions are encountered) can be described by the evolution matrix U, 
defined as follows 

b(tf) = U(tt,ti)b(ti), (10) 
b(t) = ( b b + ), (11) 

/ p -iC(tf,ti) n \ 
U(ttA) = { o e*cL))=e~ lC(t!A)a *> (12) 

u 

C(*f,*i) = \f (13) 
ti 

In particular, the phases acquired during the adiabatic stages (see Fig. [5]) are given by: 

b = l [ m<it, (2 = 1 f mdt. (14) 



B. Single passage: Landau-Zener transition 

Consider the nonadiabatic region in the vicinity of t\ 2' t = + f, uj \t'\ <C 1. Then the bias can be linearized: 

e(*i, 2 + *') w ±vt', (15) 

where 



v = Auj\cosujt h2 \ = AojJl- (^£j . (16) 

The linearized Hamiltonian 

ff(i') = -^ T T ff! (17) 

is exactly that of the LZ problem. The solution of the LZ problem in terms of the parabolic cylinder function is 
presented in detail in App e ndix A. For discu s sions of different theo retical approaches to the LZ problem, see e.g. 



presented m detail m App endix A. for discussions oi dmerent tneoretical approacnes to tne ll problem, see e.g. 
(jDelos and Thorsonl [1973) . (|Benderskii et all 120031 ). (|WittisJ . [2005T) . Here, in the main text, we also present the 



shortest solution, based on the fact that the transition under the adia batic perturbation has quas iclassical character, 
where the change of the action (given by the integral J E(t)dt) is large ( Landau and Lifshitzl . Il977h . Then the problem 



of the transition under the adiabatic perturbation is formally fully analogous to the problem of the quasiclassical over- 
barrier reflection. 

Assuming that the system is initially in the lower energy level, the transition probability from the lower to 
the upper level during the single-sweep process is described within the quasiclassical approximation as follows 
( Landau and LifshitzUl977t) . The energy levels coincide at two points in the complex plane: at t' = ±t = ±iA/v one 



finds that E+(ta) = E_(to). Accordingly, the probability P+ that the system ends up in the upper level is determined 
by the contour integral in the plane of complex time: 



exp I -21m J [E+(t') - E-{t')]dt' J . (18) 
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Calculating this integral for 2E± = ±y/ A 2 + (vt') 2 and to — iA/v, we obtain the LZ transition probability: 



Plz 



exp (— 2ir5) , 



(19) 
(20) 



Equation (Q~9|) describes the transition probability for an arbitrary value of the exponent ( Nikitinl . [l996). As the 
driving velocity v is changed from (adiabatic limit) to oo (sudden-change limit), the transition probability P + varies 
from zero to unity. In Fig. [3]the red dashed curve shows the instantaneous transition of the adiabatic-impulse model, 
with the mixing prob ability given by Eg. (1191) . The blue solid curve is calculated numerically by integrating the 
Schrodinger equation (jShevchenko et all . l2005| ) . Note that the LZ formula ([T9"l) accurately describes the final upper 
level occupation probability P + , but it does not describe the transient dynam i cs in the vicini ty of t he energy-leve l 
avoid ed crossing. This issue was studied in detail in Refs. ( Mullen et aZ.I . [l989h . ( Vitanovl . Il999lh and ( Zenesini et all 
2009). In particular, it was shown that the duration of the transition (the region between the green dotted lines in 
Fig. [3]) can be estimated as follows: 



thz ~ x v / ^max(l, VS). 



(21) 



Here we would like to note that, as argued in ( Garrawav and Vitanovl Il997f) . Eq. (|2T|) gives only the upper limit of 
the transition time, and one can think of the transition time is being shorter than the one given in Eq. (|21j) . The 
conditions of the applicability of the adi abatic-impulse model, which d escribes the LZS interferometry, require that 
the relations A/uj and/or A/uj are large ( Garrawav and Vitanovl . [l997l ): 

(22) 



A 



A 2 » u 2 . 



In particular, the multipho ton fringes were shown to be described by this model at A 
(jKrainov and Yakovlevl ll980D . 



< 
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FIG. 3 (Color online) One-passage LZ transition. The time dependence of the upper level occupation probability P+(t) is 
plotted for the following parameters: A/ A = 10, ui/A = 0.05, eo = 0. The blue solid curve presents the exact numerical 
solution, while red dashed curve corresponds to the analytical solution. The time interval between the vertical green dotted 
lines corresponds to the crossover region, which is characterized by £lz- 



The LZ formula (fT9"]) describes the upper level occupation probability P + , which is the square of the absolute value 
of the probability amplitude. However, when interference is relevant (as will be discussed in detail below), one has 
also to know the change in the relative phase between the two components of the wave function as the crossing region 
is traversed. One can keep track of this phase by describing the LZ transition using the (non-adiabatic) unitary 
evolution matrix N, which is derived in Appendix A: 



b(ti, 2 + 0) = N b(*i j3 - 0), 



(23) 



~ 7T 



(25) 



+ 5(hx5 - 1) + &rgT(l - i5), 



(26) 



where ips is the so-called Stokes phase (see e.g. ( Child! . [l974h . ( Kavanumal Il997() . ( Wubs et aZj . 12005)) and T is the 
gamma function. 



C. Double passage: Stiickelberg phase 

Consider now the double-passage process, where the avoided crossing region is passed twice at the same speed. This 
double-passage problem corresponds to one full driving cycle of the periodic driving that we are considering. The 
probability for the upper level P+ after one period is given by Eq. (|B5[) with n = 1: 

P+ = 4P LZ (l-P LZ )sin 2 $ st , (27) 
$st = C2 + <?s- (28) 

Th is result, that the probability of ex citation P + is an oscillating function of a certain phase $st , was first obtained 
in (|Landaulll932at IStuckelberelll932ft . The Stiickelberg phase $st consists of two components: the first £2 is acquired 
during the adiabatic evolution and the second tps during the non-adiabatic transiti ons. 

Stiickelberg oscillations were observed in atomic in elastic scattering cross-section s (Nikit in and Ovchinnikovalll972l ) 
and in the microwave excitation of Rydberg atoms (jBaruch and Gallagheii Il992t lYoakum et q?J . Il992h . However, in 
many cases these Sti ickelberg oscil lations average out and can be neglected, and one only needs to consider the 
averag ed probability (jNikitinl . 1 19961) 



P+ = 2P LZ (1-P LZ ). (29) 

The above expresion is the sum of two probabilities, which correspond to the system being excited at the first or at 
the second passage through the avoided-crossing region. These two components are schematically illustrated in Fig. @] 
with trajectories marked by single and double arrows, respectively. 

The quantum-mechanical interference between the different LZ transitions has the result that the net excitation 
probability after the double passage can range from (destructive interference) to 4Pl Z (1 — Plz) (constructive 
interference) according to Eq. (|27p . The latter probability is twice as large as the one wi thout interference, Eq. (1291) . 
This situation is analog ous to the one encountered in the Mach-Zehnder interferometer ( Ji et all 120031 : [Oliver et all 
l2Q05t iPezze et aU 120071) . 



D. Multiple passage 

Let us now study a multiple-passage process, where the system passes the avoided-crossing region periodically. 
Consider first the relevant time scales. The characteristic time for an LZ transition t^z can be estimated from Eq. I|21[) . 
The time between subsequent tunnelling events is of the order of half driving period T/2 = tt/lo. LZS interference 
takes place when (z) successive LZ transitions events do not overlap and (ii) the phase coherence (characterized by 
the dephasing time T2) is preserved: 

i LZ < T/2 < T 2 . (30) 

In this subsection we assume T/2 <C T2 and ignore decoherence. For definiteness, as the initial condition we consider 
the system to be in the ground state at time t\ + 0. Then the state at time t, after n full periods, is described by the 
following evolution matrices (see Fig. [2]): 



(31) 
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FIG. 4 (Color online) Double-passage transition. Adiabatic energy levels as in Fig. 2 are plotted. The lines with one and two 
arrows show the two trajectories where the transition to the upper level happens during the first and the second passages of 
the avoided level crossing. Their respective transition probabilities are given by Plz(1 — Plz) and (1 — Plz)Plz, while the 
interference is described by Eq. ([27)) . 
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(32) 



Ui.2 = exp(-i£i,2<7*). 



(33) 
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FIG. 5 (Color online) Constructive (blue dashed line) versus destructive (red continuous line) interference. Time dependence 
of the upper level occupation probability P+ with constructive (£i = 7r/2 + fc7r) and destructive (£i = kit) interference at so — 
and Plz = 0.1. The TLS is taken to be initially in the lower level. 

One can now derive expressions for the evolution matrices in terms of the parameters in the Hamiltonian (fT]), as 
explained in Appendix B. Using these matrices one can analyze the dynamics for any given initial state. Figure [5] 
illustrates two different possibilities that correspond to constructive and destructive interference between consecutive 
crossings. 1 In addition to the time dependence, one can also derive expressions for the average populations of the two 
quantum states as functions of the system parameters, since it is common to experimentally measure the steady-state 
populations of these states. In general, one needs to make approximations in order to obtain clear analytical results. 



1 These cases are also clearly illustrated in the two computer animations available online at http://dml.riken.jp/animations.php 
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Below we consider two opposite limits, depending on whether the in dividual crossings are slow (i.e., almost adiabaticj 
or fast (i.e., almost sudden). As explained in Appendix B (see also ( Ashhab et all 120071 ) ). in both cases a resonance 



condition can be derived that determines whether the system will exhibit oscillations between the two basis states. 
This condition is given by 

(1 - P LZ ) sin £+ - Plz sin £_ = 0, (34) 

where 

C+ = £i+£ 2 + 2£ s , (35) 
C- = C1-C2. (36) 

As we shall see shortly, this resonance condition results in drastically different patterns in the slow- and fast-passage 
limits. 



1. Slow-passage limit 

In the limit 6 3> 1, Plz *C 1 and the resonance condition reduces to 

C+ - Ci + C2 + 2<^s - kn, (37) 

for any integer k, and <£>s ~ — vr/2. Going further and deriving an analytic expression in terms of A, eg and A is 
complicated by the fact that the integrals that determine £1 and £2 cannot be evaluated in closed analytic form. 
Numerical calculations, however, are straightforward. The resonance condition (|37[) describes arcs around the point 
A = eo = (see Fig. |6j, where we treat A as a fixed parameter. Taking the system to be initially in the lower energy 
level, the average occupation probability of the upper level is given by (see Eq. (|B13|) in Appendix B): 

■p- = P L z(l + cos£+cos£_) , 3g . 
+ sm 2 £ + + 2P L z(l + cos£+cos£_)' 

From this expression one can see that on resonance (i.e., when sin£ + = 0) the upper level occupation probability 
is maximal, P + = 1/2. One can also see that the widths of the resonance lines are modulated by the numerator in 
Eq. (|38|) ; it tends to zero at points where 



cos C+ cos £_ = — 1 . (39) 
Figure El illustrates the interferometric pattern obtained from Eq. (|38j) . In this figure we use the parameters of 



( Sillanpaa, et all l2006t ) (see also Table fl} . Note that for any value of the energy bias offset £q , the maxima and 



minima alternate with inc reasing driving amp litude A (see the colour vari ation along the vertical dashed line in 
Fig. El), as was studied by (jShytov et aLl . l2003l ) and (|Shevchenko et aZll2005h . 



In the case of zero offset £0 = 0, one can go a little bit further with the analytic derivation. Taking £0 = we have 
Ci + C2 ~ 2A/lu and £1 — C2 = and the resonance condition (Eq. (|37p) is satisfied when 

2A/u = kn. (40) 

One can alternatively analyze the dynamics in the time domain. When £ = 0, t\ = and t 2 — k/uj, and from 
Eq. (jB5[) in the first approximation in Plz it follows that 

P+(t) = 4sin 2 CiP L z(l - Plz)^# = — - ^tf - 2 > sin2 ^> ( 41 ) 



sin 2 cos 2 Ci + Plz sin 2 Ci 



tt/uj 







where E(x) is the full elliptic integral of the second kind. 

Let us now analyze the expression for P+(t), Eq. (l4"Tj) . The prefactor of sin 2 n<p has a minimum at £1 = kir (where 
this prefactor is equal to Plz) and a maximum at £1 = -| + kir (where it is equal to 1). Thus, the former case (£1 = kn) 
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FIG. 6 (Color online) Slow-driving LZS interferometry for Auj < A 2 . The time averaged upper level occupation probability P+ 
as a function of the energy bias eo and the driving amplitude A. The graph is calculated with Eq. (|38f) for uj/ A = 0.32 < 1. The 
inclined red lines mark the region of the validity of the theory: £o < A, which means that the system experiences avoided level 
crossings. Outside of this region the excitation probability is negligibly small. The vertical dashed line shows the alteration of 
the excitation maxima and minima. 



corresponds to destructive interference and the latter case (£i = \ + kir) corresponds to constructive interference. 
These two cases are illustrated in Fig. [SJ Note that the latt er case corresponds to t he coherent destruction of 
tunn ell ing, which has been studied e xtensively in the li te rature [(|Grossmann et all 1991al lbl: iGrossmann and Hanggil . 
Il992h . (|Kavanuma and Saitol . [2008h . (|Eckel et a/.l . l2009l ). (|Ho et aZ.I . l2009h . (|Wubsl . l2010l) ]. From Eq. gl]) for the cases 
of constructive and destructive interference we obtain: 



constructive: P+{t) = sin 2 



(43) 



destructive: P+(t) = P L z sin 2 (r^nj , Ci = kw. (44) 

The dynamics (i.e., the time dependence of -P+) in the cases of constructive and destructive interference is illustrated 
in Fig . [H The resonance condition Ci = § + kir with Eq. (J42J) in the limit of strong driving (i.e., A ^> A) takes the 



form ( Shevchenko and Omelvanchou k, 2006): 

2A/uj = 7r(2fc + 1). (45) 



The positions of the resonances depend on the ratio A/u> but not on A, as was studied in detail by (|Shvtov et all . 
120031 ). Note that Eq. ((40j) predicts more resonance peaks than Eq. (|45j) . The explanation of this discrepancy is that 
at Eq — and for even values of k the resonance peaks of Eq. (|4"0")l disappear according to Eq. (|3"9")l , and we recover 
Eq. p5|) . Interestingly, if we take Eq. (|42)) in the limit of weak driving (i.e., A <§; A), we obtain: 

A = (2fc+ 1)cj. 



This relation describes the odd multiphoton resonances, as studied by (jShirlevl . Il965h and (jKrainov and Yakovlevl 



Il980h . where the energy level separation A is an odd multiple of a photon energy oj. We have therefore obtained 
the correct resonance condition, even though the weak-driving limit is outside the regime of validity of the above 
derivation, which is based on the interference of multiple LZ crossings. 

2. Fast-passage limit 

In the limit 6 <C 1, (1 — Plz) *C 1 and the resonance condition reduces to 

C- = Ci - Ca = kir. (46) 
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Unlike the slow-passage limit, one can now derive an approximate analytic expression for Ci — C2 in the large- amplitude 
limit (the reason is that the parts of the yellow and blue areas in Fig. 2 that are difficult to calculate cancel when 
taking the difference between the two areas): 



*i+2ir/(x) 



C- = Ci - C2 « 

The resonance condition is therefore given simply by 



. . . TT£q 

(eq + Asmuit) at = . 



£q = kljj. 



(47) 



(48) 



If the system is driven on resonance (with a given value k), and taking the TLS to be initially in the lower diabatic 
state (i.e., iff or tp± depending on the sign of £0)7 the system exhibits full oscillations between ip^ and tp±: 



P$(*) = £(l-cosflWt 



(fc). 



(49) 



The oscillation frequency c an be calcula t ed fr om the effective rotation angle in the full-cycle evolution matrix 
(NU2NU1), as was done in (jAshhab et aLl . l2007[ ). In the limit A > s we find that 



(k) 



A1 



cos 



A 7T 7T 

u ~ 2 " 4 



(50) 



As in the slow-passage limit, the Rabi frequency exhibits periodic behavior in the fast-passage limit as well. In 
particular, it vanishes whenever the cosine function in Eq. (|50p vanishes. The interferometric pattern that results in 
this case is shown in Fig. 7. 

Above we have used the adiabatic-impulse approach to analyze the response of the TLS to the driving field. To 
demonstrate another approach, next we consider the results of the rotating- wave approximation (RWA); for more 
details see Appendix C. 

If decoherence processes are not taken into account, the TLS is described by the Schrodinger equation. Its solution 
gives the probability of the upper diabatic state P up , which exhibits fc-photon Rabi oscillations: 



2n 



(fe)2 



1 — COS fij^ t 



A fe = AJ k ( - 



« = y/{ku,-e )* + Al 



(51) 
(52) 
(53) 



where is the Bessel function. Note that the fc-photon Rabi frequency is described by a power law at small driving, 
Ajuj <C 1; namely, at Sq = kui: 



9. 



A_ 

fc! V 2ZJ 



(54) 



Changing the system parameters (eo or cj), we pass through different fc-photon resonances. Thus, the time-averaged 
probability P up can be described by the sum of the Lorentzian-shape fc-photon resonances: 



p - Vp w - -V 



A 2 fe 



{kuj - e ) 2 + A 2 k 



(55) 



The position of the resonances is (quasi-) periodic in the parameters £q, u) , A; this corresponds to the multiphoton 
relation Eq — kuj and to the quasi-periodic character of the Bessel functions Jk{A/ui) for large value of the argument: 



•h 



-(2fc + L 



(56) 



Note that the oscillatio ns descr i bed b y Eq. ([56]) represent the counterpart to the Stiickelberg oscillations of the 
double-passage problem ( Nikitinl . 1 19961 ). 
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FIG. 7 (Color online) Fast-driving LZS interferometry for Aui S> A 2 : dependence of the time-averaged upper diabatic state 
occupation probability P up on eo/w and A/ui. The graphs are plotted using Eq. (|57|) for uj/A — 300 S> 1, loTi/(2tt) = 2.4 x 10 4 
and wT2/(2tt) = 24 (a) and w/A = 1.14 > 1, ljTi/(2tt) — ljT2/(2tt) = 6 (b). Several multiphoton resonances are shown by the 
vertical pink dotted lines at so = kuj (for k = 0, 1, 2, 5, only) modulated by Bessel functions. The vertical red double-arrow 
in (b) shows the distance between two consecutive zeros of the Bessel function. 



Note here that the results derived in the adiabatic-impulse approach and t hose of the RWA are essentially the 
same. These results can also be derived using Floquet theory ( Son et al\ , 120091 ) (see also Appendix D) or using the 



dressed-state pictur e, which represent the hybridization of the qubit and photon degrees of freedom ()Liu et Ml, [2006; 
IWilson et a/1 120071) (see also Appendix E). 

Relaxation can be described using the Bloch equations for the reduced density matrix; the influence of the environ- 
ment is taken into account using the phenomenological energy and phase relaxation times, Ti and Ti. In Appendix 
C we obtain the stationary solution of the Bloch equations, which at zero temperature results in the excitation 
probability: 

Pup — 



In the limit where the dephasing and relaxation times are large and equal to each other, i.e., T\ —Ti~^oo, the above 
expression (somewhat accidentally) coincides with the time averaged solution of the Schrodinger equat ion. Eq. (|55l). In 



Fig. FT] we plot the probability P up with Eq. (|57|) for the parameters of Refs. ( Oliver et «Ll . l2005l) and ( Izmalkov et all 
12008b (see Table H). 

We emphasize here that the distinction between the interference patterns in Figs. |B] and [7] is determined not by the 
ratio uj/A, but rather by the ratio uiy/ 'A 2 — e§/A 2 . This point is illustrated in Fig. [SI where one can see the crossover 
from the slow-passage to the fast-passage limit as the driving amplitude A is increased. Note that in order to see this 
crossover for a fixed frequency, one needs to take a small value of uj/A. 

The analysis of the positions of the multiphoton resonances allows to do spectroscopy and to obtain the parameters 
of the TLS. To calibrate the power of the generator one can use the distance between the minima or maxima in 
the Stiickelberg oscillations (one of them is shown in Fig. [7] with the vertical red double- arrow), or alternatively, the 
slope of the interference fringes (shown by the orange dashed line in Fig. [?])■ The analysis of the interference fringes 
also allows to obtain the parameters that charact erize the relaxation of the system due to the interac tion with the 
envir on ment. This was theor eti cally discussed in (IShytov et al. , 2003) and experimentally realized in [ (jOliver et all 



|2005| ). (jSillanpaa et a/.L l2006h . (jWilson et all |2007| ) . (jlzmalkov et all l2008h ] . The longitudinal T x and transverse T 2 



relaxation times can be extracted from the shapes of the resonances in Eq. (|57j) . Moreover the effective temperature 
can be estimated using Eqs. (|C15llC16| from the shape of the first resonance, where AE w |eo| ~ to which is typically 
of the order of 1 GHz and smaller for superconducting TLSs. 

In practice it can be useful to perform more detailed calculations beyond the stationary solution in Eq. (f57|) . We 
note that the relax ation times, Ti and T2, can be introduced phenomcnologically as two f itting parameters (as in 
( Oliver et al\ . \200Ei) ) or calculated within certain model for dissipating environment (as in (Sillanpaa et all l2006f) ). 



For superconducting qubits a convenient model is the spin-boson model, which models the environment as a bath of 
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FIG. 8 (Color online) Crossover from the slow-passage limit (bottom part of the figure) to the fast-passage limit (top part of 
the figure) as the driving amplitude A is increased. On the left the steady-state probability P+ of the adiabatic excited state 
is plotted as a function of bias offset Co and driving amplitude A. On the right the probability P up of the upper diabolic state 
is plotted. One can see that the resonance features are clearest in the adiabatic basis for slow passage and in the diabatic basis 
for fast passage. The ratio u/A is equal to 0.32, and there is no decoherence. Note that panel (a) differs from Fig. [6] because 
that figure was generated using Eq. (|38|l . whereas in this figure we numerically solve the Bloch equations. 



harmonic oscillators (|Makhlin et a/.l . l200ll) . (|Wendin and Shumeikol . I2007I) . (|McDermottl . l2009l) . ([Wilson et aUl20ld) : 
in this model both T\ and T 2 are parameter-dependent and determined by the strength of dissipation. However 
in many cases the relaxation can be described by effective parameter-independent times T\ and T 2 ; in this paper 
for plotting Fig. [7] we have taken the relaxation times, T\ and T2 , for s implicity being constant. Also, the Bloch 
equations (IC13I) can be solved exactly as, e.g., in (Sh evchenko et aLl . l2008l ). which allows one to relax the assumption 
A <C u>; the positions of the multiphoton resonances are then defined by the relation: AE ps A 2 + e§ rj ku. 
Note that the inhomogene ity of the field due to low-frequency noise cannot be directly described within the Bloch 
equations ( Abragaml . 1 1 96 if) ; it results in the broadening of the Lorentzian-shaped resonances and can be described as 
the additional term in the expression for the width of the resonances, as in Ref. (|01iver et aLl . l2005l) . 



E. Decoherence 

When studying the effect of decoherence on the dynamics of a quantum system, it is perhaps most common to start 
from the case of zero decoherence and slowly increase the decoherence rates. In the problem of LZS interferometry, 
we find it simpler to think about the effect of coherence rather than the effect of decoherence. We study the effect 
of coherence by co nsidering the steady-state populations of the two quantum states of the TLS (as was done in 
(|Berns et all [2006) ). We note here that in the slow-passage limit we analyze the populations of the adiabatic-basis 
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states whereas in the fast-passage limit we analyze the populations of the diabatic-basis states. 
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FIG. 9 (Color online) Same as in Fig. [6] (i.e. LZS interferometry with low- frequency driving), but including the effects of 
decoherence. The time averaged upper level occupation probability P+ was obtained numerically from the Bloch equations 
with the Hamiltonian JT]). The dephasing time T2 is given by ujT2/(2n) — 0.1 in (a), 1 in (b), 5 in (c) and T2 = 2Ti in (d). The 
relaxation time is given by u)T\/(2n) = 10. 



Starting from the limit of strong decoherence, no interference between the different LZ transitions occurs. The 
result is then simple. If the driving amplitude is large enough to take the TLS through the avoided crossing, each 
crossing will act as a small "kick" to the populations of the two states, and these kicks will eventually add up such 
that the TLS reaches a steady state with equal population of the two energy levels (the ground state can have a 
higher population if relaxation is substantial). If the driving amplitude A is smaller than the bias offset \eq\, no LZ 
transitions will occur, and the system will remain in the lower level. This situation is depicted in Figs.^a) andQjIJa). 

We now increase the coherence time such that the two transitions in a single driving period are separated by a 
time that is smaller than the coherence time. The transitions induced by different driving cycles are separated by 
more than the coherence time, and they therefore act as independent kicks. As in the case of strong decoherence, a 
sufficiently large amplitude is required to drive transitions between the two quantum states. However, if the phase 
accumulated between the two interfering transitions corresponds to destructive interference, the net effect of a full 
driving period will be no mixing at all. As a result, we now find lines in the A-Eq plane where the TLS is driven back 
and forth across the avoided crossing but the system remains in the lower level (in the fast-passage limit, the TLS 
remains in the diabatic state that corresponds to the lower level for the given value of Eq). This situation is depicted 
in Figs.^b) andHUIb). 

We now increase the coherence time further, such that interference between transitions from a few successive driving 
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cycle s occurs. [Relate d discussions of the intracycle and intercycle interference were given in (lAshhab et all l2007h 
and (|Arbo et aUl2010t) ]. We now start to see the resonance lines forming, or rather the mixing between the two states 
being suppressed whenever the resonance condition is not satisfied (Figs. EJc) andHUf c)). Finally, in Figs. (9jd) and 
HUH ) we take the case of negligible decoherence, and we find that the resonance lines are now well defined. Deviation 
from the resonance condition leads to the absence of any substantial mixing between the two states. 
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FIG. 10 (Color online) Same as in Fig. [7] (i.e. LZS interferometry with high-frequency driving), but including the effects of 
decoherence. The time-averaged upper diabatic state occupation probability P up is obtained numerically by solving the Bloch 
equations with the Hamiltonian fl}. The dephasing time T2 is given by uiT2/(2n) = 0.1 in (a), 0.5 in (b), 1 in (c) and T2 = 2Ti 
in (d). The relaxation time is given by ujTi/(2-k) — 10 . 



F. "Breakdown" of the adiabatic theorem 

The adiabatic theorem i s used in various subfields of physics. The basic idea of the theorem is the following 
( Landau and Lifshitd . Il977h : if a physical system starts in one of the eigenstates of the system's Hamiltonian, the 
Hamiltonian is varied slowly enough and there are no degeneracies between the eigenstates during the variation of the 
Hamiltonian, then the system will end up in an eigenstate of the Hamiltonian whose order (starting from the ground 
state counting upwards) is the same as that of the initial state. 

As the term "theorem" suggests, the above statement is a rigorously provable statement. The tricky point is in 
deciding how slow the Hamiltonian variation must be in order to ensure adiabaticity. The most common form of the 
adiabatic condition is obtained by inspecting the LZ formula (Eq. (TH)1) ): 6 ^> 1, and slightly generalized forms of this 
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condition. If the condition 5 3> 1 is satisfied in the LZ problem, the system will almost certainly stay in the ground 
state, assuming it started in the ground state. The problem is that the above condition was derived for a specific 
problem, namely the simple LZ problem, and not from general considerations that apply to any quantum system with 
any f or m of Hamiltonian vari ati on. This issue has been t he subject of a n umber of recent papers [ ([Marzlin and Sanders! . 
12004 ), (jSarandv et all |2004) , (|Tong et all 120051 12007() . (jAminl . l2009l) ]. For example, if a many-level system starts 
in an eigenstate of the Hamiltonian, and then during the variation of the Hamiltonian an enormously large number 
of avoided level crossings are encountered (in w h at can colloqu ia lly be called a "spag hetti" energy- level structure) , 
as e.g. in [ (jPokrovskv and Sinitsvnl . 120021 12004) . (|Shvtovl . |2004) . (jGaranin et all 120081 )] . knowing that each crossing 
depletes the population of the initial eigenstate by a small amount does not guarantee that the total depletion of the 
state from all the crossings will also be small. 

Another example where the above simple condition does not guarantee adiabaticity over the full time evolution 
occurs in the simple LZS interferometry setup considered in this review. If the system is repeatedly driven back 
and forth across the avoided-crossing point and the phase accumulated between successive crossings is designed to 
produce constructive interference, the small transition amplitudes of the different crossings will add up to produce a 
large transition probability. For example, if we take S 3> 1, sq — and A/uj = 7r/2 + A:7r (i.e., we choose the parameters 
to satisfy the resonance condition, Eq. (I43[l ). we know that it is only a matter of time until the ground-state population 
is transferred completely to the excited state and vice versa. Finally, it should be emphasized that one does not need 
to satisfy any resonance condition of the kind given in Eq. (1371) . An arbitrary- looking driving signal that produces 
suitable phase accumulation between successive crossings will also eventually result in full population transfer between 
the ground and excited states. Even random phases will result in diffusion-like dynamics that ultimately result in 
large mixing between eigenstates. 

When discussing the adiabatic theorem, it is worth digressing for a moment to mention an ide a that has received 
considerable attention recently, namely adiabatic quantum computing (AQC) (jFarhi et all l200ll) . The basic idea of 
AQC is as follows: one subjects a physical system to a simple initial Hamiltonian, such that it is guaranteed that the 
system will relax to its ground state. One then adiabatically changes the Hamiltonian to reach a final Hamiltonian 
whose ground state encodes the answer a given problem (the specifics of the problem are encoded in the parameters 
of the final Hamiltonian). Measuring the final state then provides the answer to the problem. During the evolution 
from the initial to the final Hamiltonian, one can expect that avoided crossing structures will be encountered, and the 
size and number of these avoided crossings will determine the requirements on the sweep rate. The smaller the energy 
gaps and the larger the number of crossings, the more slowly one has to sweep the parameters of the Hamiltonian. 
The traversal of the avoided crossings is governed by LZ physics, and non-adiab atic LZ transitions pl ay a destructive 
role and should be avoided; this was studied for a superconducting qubit e.g. in ( Graicar et aU 120051) . It is currently 
unknown how the running time scales with system size for different problems. The answer to this question will 
determine whether AQC provides a computational advantage ov er classical computer s or not. In a study that comes a 
step closer to the subject of the present article, it was argued in ( Ashhab et aU 120061) that the interaction between an 
AQC system and its surrounding environment could lead to a situation where the environment effectively drives the 
system back and forth across the avoided crossings, thus reducing the probability of remaining in the ground state. 



G. Reversal of stimulated emission 

It is well known that when an atom is placed in a resonant cavity and the atom is initialized in its excited state, 
it emits a photon into the cavity with a rate that is proportional to the number of photons in the cavity (more 
accurately, the emission rate is proportional to n + 1, where n is the number of photons in the cavity). Recently, 
this situation was realized in an experi ment that uses a sup erconducting charge qubit whose energy-level structure is 
essentially the one that we study here ( Astafiev et all [20071) . The qubit is biased such that it experiences an effective 



inverted relaxation taking it from the ground to the excited state. The qubit therefore constantly emits photons into 
the resonant cavity , and a lasing s tate i s realized. 



As discussed in (jAshhab et all l2009i) , for small photon numbers the emission rate is proportional to the number 
of photons in the cavity, which can be seen as being a result of the fact that the Rabi frequency is proportional to 
the square-root of the photon number. It is interesting to consider what happens if the photon number is so large 
that the cavity field effectively drives the qubit back and forth across the quasi-crossing region. In this case the Rabi 
frequency is no longer proportional to the driving amplitude (which is proportional to the square-root of the photon 
number), but it follows the Bessel- function dependence that is characteristic of LZS interferometry (see Appendix 
C). One therefore has the situation where the emission rate is no longer proportional to the photon number in the 
cavity n, but rather it is proportional to J^{ot\/n), where a is a conversion coefficient. Stimulated emission therefore 
exhibits oscillatory behavior, where the emission rate has regions where it increases with increasing photon number 
and regions where it decreases with increasing photon number. 
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One consequence of the quasi-periodicity in the Bessel function is that, with a suitable choice of parameters, there 
can exist several stable steady-state photon numbers in the cavity. These solutions can be obtained by comparing the 
emission rate into the cavity (which follows the Bessel-function dependence) and the loss rate out of the cavity (which 
is proportional to the phot on number). Thi s situat ion concerning the existence of multiple steady-state solutions in 
this system was studied in (jRodrigues et all l2007alfbh . 



III. RECENT EXPERIMENTS WITH STRONGLY DRIVEN SUPERCONDUCTING QUANTUM CIRCUITS 
A. Multiphoton transitions in superconducting qubits 



If a quantum system is strongly driven, it can be excited by a multiphoton process. One way to observe the 
multiphoton resonanc es is to con s ider t he frequency dependenc e, which is con v enient for atoms exp osed to a periodic 
electromagnetic field [ (jGallagherl Il994r i. (jFregenal et aZ.I . l2004h . (|Fe>rreL \200A) . (jMaeda et aLl . l2006h ]. Another way is 
to tune the energy difference at fixed driving frequency. In superconducting qubits the energy difference is tuned by 
an external source. The observation of the resonances allows for multiphoton spectroscopy. Furthermore, the strong 
driving provides a tool to control the energy level population, with the driving frequency being a fraction of the energy 
level difference. This is important for operating superconducting qubits, where the energy differen ce is of the order 
of 0.1 to 10 GHz. The multiphoton resonances have Lorentzian shapes ( Tornes and Stroudl 120081 ). as described by 
Eqs. (I55p or (|57[) . which provide a way to estimate the relaxation rates in the system. 

Multiphoton resonances were observed in sever a l supe rconducting qubits with a small number of photons. They 
were observed in charge qubits ( Nakamura et all l200lf ). where the energy- level difference is controlled by a gate 
voltage, as well as in the c harge-phase qub i ts, where th e energy difference is controlled by both a gate voltage and 
an ext ernal magnetic flux ([Shnvrkov et all 120061 l2009f) . Up to five-photon transitions were reported for the phase 
qubit (|Wallraff et al\ . \200S) . where the energy difference is tuned by the bias current. Multiphoton transitions were 
also observed in a flux qubit ( Saito et all 12004 ) controlled by an external magnetic flux. Further improvement of the 
coherence characteristics of the superconducting qubits allowed to observe resonant excitations with more photons, 
at higher driving amplitudes. These results are presented below. 



B. LZS interferometry in superconducting qubits 



The observation of LZ tunnelling in supe rconductin g qubits was first reported by (Izmalkov et al\ , \2004) and then by 
( Ithier et aU fepOS: Jo hansson et aZ.I . I2009T ) (see also ( Ankerhold and Grabertl . 12003 ) for the theoretical background). 
Recently LZS interferometry was demonstrated by se v eral groups on differen t superconduct i ng qu b its with comple- 
ment ar y measuremen t techn iques [( Oliver et aLl . l2005l ). ( Sillanpaa et aLl . l2006l ). ( Wilson et aZ.I . l2007h . ( Izmalkov et all 
I2008T) . (jLaHave et all . 12009)]. Below we briefly describe these works and present the major results in Fig. [TT] There, 
the left column shows the circuits used in the experimental setups and the right columns shows their corresponding 
experimental LZS interferometry. The latter shows a measurable quantity versus driving power (amplitude) and en- 
ergy bias (which de fines the energy l evel distan ce). Also recentl y the LZS interferometry was demonstrated with the 
rf SQUID qubit in (|Sun et all [2009) [see also (|Sun et all 120061 ) for experimental details] . One important difference 
from the previous works is that the system is essentially a multi-level system with several different avoided crossings 
involved (see also next subsection). This allows for additional interesting and useful phenome na such as contr ollable 
population inversion, as convincingly demonstrated both experimentally and theoretically in ( Sun a/] . l2009h . 

The parameters used in different experiments studying LZS interferometry with driven superconducting qubits are 
given in Table |TJ Remarkable is the fact that in all of these experiments the relation of the characteristic time of the 
LZ transition (Eq. (f2"TjO estimated with maximal driving amplitude ^4 max 5 t LZ ~ l/-\A4 max uj, to the driving period 
T was almost identical: tuzjT ~ 0.05. Note that the exact value of t^z/T does not affect the analysis or results, 
provided that this ratio is much s maller than one (see Eq. (|3"0"|) ). 

(1) In Ref. ( Oliver et all l2005h a superconducting flux qubit was subjected to a strong microwave drive and LZS 
multiple interference fringes were observed. The flux qubit was a niobium superconducting ring with three Josephson 
junctions. The circuit and results are shown in Fig. [TO a) to the left and to the right respectively. The qubit 
was biased with both dc and microwave magnetic fluxes. The level population was probed with a dc SQUID. The 
switching probability of the SQUID to the resistive state measures the state of the qubit in the diabatic basis {</?-)-, <Pi}- 
These measurements showed resonant peaks and dips when the energy of k photons matched the qubit 's energy- level 
difference. The latter was controlled with the dc flux (flux detuning) and multiphoton resonan ces with k up to 2 
were observed. Later the observation of multiphoton resonances of order up to 45 was reported (|Berns et a t, 2006). 
The interference fringes shown in Fig.[TTJare described by the theory presented in Sec. II. D. 2; cf. also Fig.JTJa). The 



20 



TABLE I Parameters used in different experiments studying LZS interferometry: tunnelling amplitude A, maximal driving 
amplitude J 4 max ] and driving frequency u) in the units GHzx27r, minimal adiabaticity parameter 5 mm = A 2 /(4uM max ), and 
maximal LZ probability P^ ax = exp(-27r(5 min ). 





A 


^raax 




gmin 


r>max 
■TLZ 


(Oliver et al, 2005) 


0.004 


24 


1.2 


w- 7 


1 


fSillanuaa, et al, 2006) 


12.5 


95 


4 


0.1 


0.5 


(Wilson et al, 2007) 


2.6 


62 


7 


0.004 


0.98 


(Izmalkov et al, 2008) 


3.5 


40 


4 


0.02 


0.9 



fringes exhibit the Bessel- f unction dependence , according to Eqs. ([52| l55" j) . and the steplike pattern in Fig. flTTa) was 
called a "Bessel staircase" ( Oliver et ali [2005h ; the inclined white lines in Fig . fTTTa) show the ma xi ma for Jl (A/uS) 
which were numbered with the Roman n umbers I to VI. For more details see (jRudner et ali [2008T ) . ( Bvlander et ali 
l2009t ) and (lOliyer and Valenzuelal. 120091) . 

(2) In Ref. (Sillanpa a et ali 120061 ) an interferometer-type charge qubit was studied. The qubit consisted of two 
nearby Josephson junctions embedded in a loop (Fig. Illf b)). The charge on the island between the junctions (the 
so-called Cooper-pair box) was controlled via the applied gate voltage V g — en g /C g , n g (t) = n g o + Sn r { sin(u; r fi). The 
effective Josephson energy of the two junctions was tuned using the applied magnetic flux $ through the loop. The 
energy levels of the qubit were controlled by the flux $ and by the dc component (n g0 ) of the voltage, while the energy- 
level populations were controlled by the ac component of the voltage. The qubit was coupled in pa rallel to an LC- 
circuit. The effective capacitance C c g of the qubit modifies the resonant frequency of the LC-circuit ( Sillanpaa et ali 
2005). This provides a method for the continuous-time monitoring of the qubit 's effective capacitance, which is 
related to the occupations of the qubit's energy eigenstates. The resulting LZS interference pattern is shown to the 
right in Fig. [TTTb) ; thi s is described by the t heory presented in Sec. II.D.l; cf. also Fig. [6j For more details see 



also ( Paila et ali 120091: ISillanpaa, et al , 2007). Also v ery recently LZS interferometry with an interferometer- type 



charge qubit wa s demo nstrated in ( LaHave et ali I20091) . The parameters of the qubit were close to the ones in Ref. 
(Sillanpa a et ali [2006); consequently, LZS interference fringes were similar to the ones shown in Fig. ITTT b). The 
important difference is in the measuring device : the qubit was coupled to a nanoelectromechanical resonator, as 
studied theoretica lly in ((irish and Schwab]. | 2003l) . and the nanoresonator frequency shift was measured. 



(3) References (jWilson et al. 



20071 120101) studied a single Cooper-pair box (SCB), which was composed of an 
aluminium island connected to a reservoir via a Josephson junction (Fig. 111( c)). The state of the qubit was controlled 
by the gate voltage V g (n g = C g V g /2e) and by the magnetic field in the dc SQUID loop connected to the island, which 
tunes the Josephson energy Ej. The q ubit was driven through the gate with the microwave amplitude A^. The qubit 
state was probed with an rf oscillator ()Dutv et al .], 120051 : iPersson et al\ , [2010) . The color image in Fig. Illf c) presents 
the phase of the rf reflection coefficient T. The ground and excited states contributed phase shifts with opposite 
signs. The result s were interp r eted i n terms of the dressed states of the qubit coupled to the driving microwave field 
dLiu et a/.l .l2 006D. (lGreenberkl2007h . 

(4) In Ref. (Izmalkov et aLl . l2008[ ) a superconducting flux qubit was driven by an external magnetic flux, analogously 
to Ref. (jOliver et ali 120051 ) . The qubit was weakly coupled to a classical resonant tank circuit, Fig. ITlTd). The tank 
circuit was biased with an rf current /bias and the measurable quantity was the phase shift 6 between the voltage 
V and the bias current in the tank circuit. The inductive coup ling of the qubit to th e tank circuit r esulte d in a 
phase shift determined by the effective inductance of the qubit ( Greenberg et ali |2002| ; lll'ichev et al. 1. 120091) . The 
resonant excitation of the flux qubit changes the direction of the current in the qubit's loop, which resonantly changes 
the effective inductance of the qubi t. The respective r esona nces were observed in the phase shift dependence as a 
sequence of the ridges and troughs ( Shevchenko et ali I2008T ) . The LZS interference pattern is shown in Fig. ITTTd) 
in the dependence of the tank phase shift on the dc (<i>b) and ac (^ac) components of the magnetic flux; the 
interference fringes are described by the theory presented in Sec. II: see Fig. [T^b) . 

Thus, strongly driven superconducting qubits were studied by different, complementary methods; their resonant 
characteristics were shown to be periodic in the bias Eq and driving amplitude A due to Stiickelberg oscillations. 
The different experimental measurements allowed to obtain the following parameters: (i) those which characterize 
the qubit and enter in the pseudo-spin Hamiltonian, (ii) the driving microwave amplitude as felt by the qubit (which 
is the calibration of a generator signal), (iii) the relaxation and dephasing rates, which describe the impact of the 
environment on the qubit dynamics (via the Bloch equations). 
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C. Multiphoton resonances in multi-level systems 



Most qubits are formed by the lowest two levels in a multi-level structure. Moreover, a qubit can be coupled 
to another quantum system, in particular, to another qubit, so that the overall structure is also multi-level. High- 
amplitude driving allows not only to reach the upper level of a TLS (qubit) with the absorption of several photons, 
but also to involve more of the qudiVs energy levels in the process. 

Multiphoton resonant excitations were studied in superconducting multi-level systems. One-, two-, and three- 
photon transitions were driven in the flux qubit not only between the two lowest levels (which form the qubit), 
but rather between five levels in ( Yu et al In this case multiphoton spectroscopy allowed to obtain the 

system paramet ers and to vi s ualize the energ y-level diagram. The l ow- frequency and ex tremely strong-driving regime 
was studied in (jBerns et all 120081 ) [see also (|Wilhelml . l2008h and (|Wen and Yul . l2009h ]. Berns et al. demonstrated 
that their method, which they called amplitude spectroscopy, allows probing higher levels that are three orders of 
magnitude higher in energy than the driving frequency. Recently interference fringes associated with LZ transitions 
at two nearby avoided crossings were also obse rved in a multi-lev el structur e of the rf-SQUI D qubit driven by large- 
amplitude and high-frequency microwave field (jWang et «Ll . l2O10t) [see a lso dSun et al , 2009f) for more details]. Also, 
multiphoton multi-level Rabi oscillations were studied in a phase qubit ( Dutta et all 20081 ). Multiphoton transitions 
in the multilevel system formed by the dr essed states of a qubit coupled to a quantum resonator were demonstrated 
in ( Bishop et all 120091 : iFink et all , I2009T ) . Another situation where a multi-level structure appears, is a system of 
co upled qubits. Direct and ladder-type multiphoton excitations were recently demonstrated in a two-flux-qubit system 
in (|H'ichev et a/.l . l2010fh 



IV. CONCLUSIONS 



Strong driving of a TLS results in a periodic dependence of its response on the bias, driving frequency and ampli- 
tude. This dependence was recently observed in su perconducting qubi t systems by several gro ups using complemen- 
tary measurement tec hn iques, both in flux qubits (|01iver et aZ.I . l2005h . (jlzmalkov et all 120081 ) and in charge qubits 
(Silla npaa et al. . 120061 ). ( Wilson et aZ.I . [2007h . Very recently LZS interferometry was also realized in a rf SQUID qubit 



Sun el all 12009). No related effect has been reported in experiments using phase qubits. This is because it is difficult 



to produce avoided energy-level structures in phase qubits without involving higher levels. 

LZS interferometry allows, first, to control the qubit state, and, second, to obtain the system parameters. In 
particular, from the positions of the multiphoton resonances, the qubit parameters can be obtained (spectroscopy); 
the Stiickelberg oscillations can be used for the definition of the driving amplitude (calibration of power); and the 
width and shape of the interference fringes can be used to estim ate the relaxation ra tes. These properties make LZS 
interferometry a very powerful tool for studying TLSs (see also ( LaHave et all [2009) . where LZS interferometry was 



realized recently for a superconducting qubit coupled to a nanoelectromechanical resonator) . 

We have reviewed recent experimental results on LZS interferometry, and we have presented the theoretical de- 
scription of this problem. In particular, we have demonstrated the applicability of the adiabatic-impulse model for 
both slow and fast driving regimes. This model is a powerful tool in the study of non-adiabatic transitions at avoided 
crossings and the interference between these transitions. 
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Appendix A: Solution of the Landau-Zener problem 
1. Adiabatic wave function 

In this Appendix we consider the driven two-level system in the adiabatic representation, assuming the evolution 
to be adiabatic far from the avoided level crossings, with non-adiabatic transitions occurring around these points. 
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The wave function in the adiabatic approximation takes the form 

rj> = ^ °* V ai, 2 = A 1)2 exp ^-t J E(t)dtj . (Al) 
Then the Shrodinger equation if^f? = £fy gives the instantaneous eigenvalues, Eq. and eigenvectors: 

V>± = <p±e*« + *\ C = ~ / n(t)dt, (A2a) 
/3 T ^ a i /n±e 



^-U* J (A2b) 

Here the term 7r/4 arises w hen the next term in the quasiclassical (adiabatic) approximation is taken into account 
( Delone and Krainovl Il985l) . 

Let us consider the original LZ problem described by the Hamiltonian ([1]) with e = vt {see also Eq. (ITTl) ; we consider 
t instead of t' in this Appendix}: 

H(t) = ~ ^ ~ (A3) 

Consider now the adiabatic wave function in the region far from the avoided level crossing, where |e(t)| = v \ t\ 3> A. 
Note that this latter condition is consistent with the condition uj \t\ <C 1 used for the linearization of the harmonic- 
driving problem, Eq. (|15|) . at large enough driving amplitudes. Then we have: 

» [ 1 ], ¥>-«[?] if*a<0, (A4a) 



o / Vi 













i 




1 







if i, > 0, (A4b) 



C(±0=^ /n(t)*«± 



(A5) 



The subscript "a" indicates that we now consider the time within the asymptotic region, where \e(t)\ ^> A, and the 
acquired phase was split into time-dependent and independent parts: 



*W = y + ^ln(V2z a ), 



4>A 



~6(lnS-l). 



(A6) 
(A7) 



2. Non-adiabatic transition (evolution matrix for the Landau-Zener transition) 

Let us now consider th e non-adiaba tic transition between the energy levels in the vicinity of the point of avoided 
level crossings, following ([Zenerl . 119321 ). Then the Schrodinger equation itp = Htp takes the form: 



This can be written in the form of the second-order Weber equations: 



(A8) 



^^ + (2STi + z 2 )a lt2 = 0, (A9) 



V -t. (AlO) 
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The solutions are combinations of parabolic cylinder functions (jGradshtevn and Rvzhikl . Il994h : 

ai = ^A±Z?_i_ 45 (±V2e l7r/4 z), (Alia) 
± 

a 2 = ^S ± L»_ 45 (±V2e l7r / 4 z). (Allb) 
± 

We find the relation between the coefficients A± and B± from the first-order equation (by inserting the solutions for 
ai t 2 in the first equation of the system (|A8p ) making use of the recurrence formula: 

B± = I f—j=—A±. (A12) 

Now with the asymptotes of the parabolic cylinder functions (|Gradshtevn and Rvzhikl . Il994f ). we find far from 
the transition point (where z = 0), at z = ±z a with z a 3> 1: 

ai (-z a )«A + S ie **H (Al3a) 

a 2 (-z a )n(-e-i 5 A + + e^A^)E 2 e^H (A13b) 

ai (z a )«v4_S ie l *( 2 ^, (Al3c) 

a 2 (z a )n(-ei s A + + e -^A_)S 2 e~ J * ( H (A13d) 

where 

E 2 , -Lexp(-i|-^). (A14b) 
Let us match this asymptotic solution with the adiabatic one, ip(t) = b±<p±, from the previous subsection: 

at t a < 0, (Al5a) 

at t a > 0. (A15b) 

Here the time dependence in the l.h.s. and in the r.h.s. is the same, see respectively Eqs. (|A5[) and (IA13|) . This allows 
to describe the evolution from t — —t a to t a as a sequence of adiabatic evolutions from t = — t a to t = —0 and from 
t = +0 to t — £ a , with the nonadiabatic transition in between (described by the time- independent matrix N), namely: 

b(t») = U(t a , +0)NU(-0, -t a ) b(-t a ) = e-^^Ne-^M-Q- (A16) 

Then collecting this with Eqs. (|A15|) and (|A13j) and equating multipliers with A + and A-, we obtain the expression 
for the matrix N. This matrix describes the transition between the two adiabatic states with the probability 
given by Eq. (fl9|) and also includes the phase jump 1^3. We note that the matrix N at t = t 2 is the same as at t = t\ 
for the nonadiabatic transitions. 

The argument of the asymptotes of the gamma function T are ( Gradshteyn and Rvzhikl . 1 1994T ) : 

argr(1 -^ )ss { m a uT^\ (A17) 

[ -4 - S{lnS - 1), 5 > 1, 

where C ss 0.58 is the Euler constant. Thus, the phase jump fs{S) is a monotonous function, which changes from 
in the adiabatic limit (S 3> 1) to 7r/4 in the diabatic (fast driving) limit (S <?C 1). 

It is worth making a comment on bases here. In this review we have used the energy eigenbasis for describing the 
dynamics of the TLS. Furthermore, we have followed the reasonable convention where the ground state far to the 
right of the avoided crossing region coincides with the excited state far to the left of the avoided crossing region, and 
vice versa. In principle, we could have defined these pairs of states such that they differ by some phase factors. If 
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we follow such a convention, there would be some asymmetry between the different matrices (denoted by N) that 
describe the different LZ transitions. However, since such a convention seems rather unnatural, we do not comment 
on it further. A different convention from ours that appears more frequently in the literature is the use of the diabatic 
states for describing the dynamics. In this case the matrices that describe the evolution from time — t a to time i a for 
a single passage depend on the signs that appear in the Hamiltonian, Eq. (|1 T|) . The part describing the dynamical 
phase accumulated during the adiabatic parts of the evolution is straightforward. The LZ transition matrices are 
less obvious, and we give the rules for how to determine them explicitly here. For the Hamiltonian in Eq. (|17p with 
negative signs in front of both terms in the Hamiltonian, the LZ transition matrix takes the form: 



(A18) 



If the sign in front of the second term in the Hamiltonian is reversed, the off-diagonal matrix elements in Eq. (|A18|) 
are reversed. If the sign in front of the first term in the Hamiltonian is reversed, only the minus sign at the beginning 
of the bottom-left matrix element is moved to the top-right matrix element. Combining the two rules above, if both 
signs in the Hamiltonian are reversed, only the minus sign in front of the Stokes phase is moved from the bottom-left 
matrix element to the top- right matrix element. Note that in the diabatic basis there will always be two relevant 
signs in front of the second term in the Hamiltonian, depending on whether the avoided crossing is traversed from 
left to right or from right to left. As a result, there will be two different matrices describing the different crossings. 

Appendix B: Evolution of a periodically driven two-level system 

The time evolution of a periodically driven TLS is described by Eqs. ([3T]) and (|32|) . Consider first the matrix for a 
one-period evolution: 

NU 2 NU 1 = ^ ° "J ) , (Bla) 

a = (I- P L zK iC+ - PLze- lC - , (B2a) 
7 = V^lz(1 - P L z)e^ s (e-<+ +e-*-), (B2b) 



C+ = (i+b + ^s, (B3a) 
C- = C1-C2. (B3b) 

Af ter diagonalizing the matrix, it is straightforward to write down its n-th power, e.g. ( Bvchkov and Dvkhnd . 
Il970l ). This unitary evolution matrix (NU2NUi) n has the form 

{NU 2 NU 1 ) n = I " n ~ U * 21 ) , (B3c) 

with elements: 

sin n(j> . 

Mil = cosncp + i(lma) — : , (B4a) 

sin0 

sinnd) . , , 

U21 = 7^-S B4b) 

cos</> = Rea. (B4c) 
Then with Eq. (|3I[) for (t — 2irn/ui) € (£1,^2)1 we obtain the upper-level occupation probability, P + = |6 + | 2 , 

□(I)/ \ 1 ,2 sin 2 n0 

P+(n) = 7 . 2 7 ■ ( B5 ) 
sin 
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And for (t - 2nn/uj) G (t 2 ,h + 2tt/uj) with Eq. j32J we obtain: 

„ni'i / \ „ sin 2 nd> _ sa\2n<h „ ,„ , 

Pr ) {n) = 2Q 1 v+Q 9 ——Z.+p LzCaB 2n4>, B6 

sin <p sin 

Qi = Plz [Plz sin 2 C- + (1 - ftz)(l + cosC+ cosC-)] , (B7a) 

g 2 = 2P L z(l - Plz) cos(Ci + fa) cos(C 2 + fa). (B8) 
Now the time-averaging corresponds to averaging over n^> 1 periods and, from Eqs. (|B5|) and (|B6|) . we obtain: 



pd) = hi 2 = i bf fBg s 

+ 2sin 2 2| 7 | 2 + (i ma )2' 



p (ii) = (B1Q) 
sin 4> 

From Eq. (IB9[) it follows that the upper-level occupation probability is maximal at Ima = 0. This results in the 
general resonance condition: 

(l-P LZ )sinC+-P L zsinC- =0. (Bll) 

In particular, in the slow and fast limits the resonance condition takes the form: 

Ci + C2 = kir for 2tt8 > 1, (B12a) 
Ci - C2 = kir for 2nS < 1. (B12b) 

The resonance condition (|Blip was derived in (|Ashhab et all |20Q7[) using an alternative, geometric visualization 
of the dynamics. There, the evolution matrix for one full driving period, i.e. (NU2NU1), is decomposed into two 
rotations, a rotation around an axis in the xy plane followed by a rotation around the z-axis. The resonance condition 
is then identified as the requirement that the z-axis rotation does not affect the dynamics, i.e. the z-axis rotation 
is a rotation by a multiple of 2n. Straightforward algebra then gives Eq. (jBllI) . A related alternative derivation of 
the resonance condition relies on visualizing the evolution matrix for one full driving period as a single rotation: on 
resonance the rotation axis lies in the xy-plane 

First, let us consider the case of the slow-passage limit, where <5> 1 and Plz *C 1. We now neglect the difference 

between pP and Pj 11 "' and, from Eq. (|B10j) to the lowest order approximation in Plz, we obtain: 



PT= 2 JW + ^C+oobC-) . (B13) 

sin C+ + 2Plz (1 + cos C+ cos ) 

For small offset So/A <C 1 the expressions fo r C+ can be simplifi ed: Ci + C2 ~ 2A/ui, (1 — £2 ~ neo/co. This makes 
Eq. (1B13[) analogous to the result obtained in ( Shvtov et a t, 2003) (their Eq. (20) at u = w —> 0). But for arbitrary 
(not small) offset £0 one has to make use of Eq. (IB13I) with (± given by Eqs. (|B3[) and (Tl4")) ; the result, such as Fig. [5] 
is significantly dependent on this correction. 

One can also obtain an idea about the form of the resonance condition in the zeroth in A /sq approximation: 



1 



Ci « - / e(t)dt 



_ 7T£ 


£o 




cos 


UJ 


UJ 








cos — 


UJ 


A 



A uj 



£0 



2 

(B14a) 
(B14b) 



Equation (|B12a[) then gives 



£0 + 2 VgHfi _ 2 £^ cos _i 50 = k _ (B15) 

UJ 7TUJ 7TUJ A 
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We emphasize, however, that the above equation should be seen as a limiting case for small A that demonstrates the 
complicated form of the resonance condition. 

Now consider the case of fast passage, where (1 — Plz) ~ 27r5 <C 1. The fast limit (8 <C 1) assumes a small value 
of A. In this limit there is a large probability (Plz ~ 1) for transitions between the adiabatic states in a single- 
passage; while the transition probability between diabatic states is small, (1 — Plz) *C 1. Hence, we will consider the 
time-averaged probability of the upper diabatic state P up . Note that upper diabatic state is ip^ for £o < and ip\, for 
£o > 0. Let us consider the latter case, as shown in Fig. [TJ Then we can relate P up to pf^ and pj 77 ' 1 as follows: 

Pu V ~P { + ] for (t-27rn/w) e (ti,t 2 ), (B16) 



P up «l - P{ 7/) for (t - 2im/u) G (h,h + 2n/ui). (B17) 



Then we note that 1 — P+' ~ P+ and obtain: 



1 4cos 2 (C 2 - tt/4) 

2 sin 2 C- +4cos 2 (C 2 - tt /4) ' 



^up^o - "T7 — . - ---9//- 777T' ( m °) 



where £ 2 is given by Eq. (|B14b[) . On resonance £_ = fc7r and P up = 1/2, then with ps tteq/uj for small offset 
*C 1) w e obtain the resonance frequency a/" = So/k. Consider now the vicinity of the fc-th resonance, 
introducing a small frequency detuning as follows 

8uj = kuj- e . (B19) 

Then to first approximation in So/A, we obtain: 

< = <B2 ° a) 

A fe = AW^-cos f - - ^(2fc + 1) ) . (B20b) 



ttA \ui 4 

The total probability P up is obtained as the sum of the partial contributions P up . This result, obtained in the 
adiabatic- impulse model, coincides with the one obtained within the RWA, Eqs. (|55H56I) . 



Appendix C: Rotating-wave approximation 
1. Hamiltonian in the RWA-form 



Following (jPegg and Series! . 1 1973ft . (|Lopez- Castillo et aU [1992ft . (|01iver et all |2005[ ) and (|Ashhab et all 12007ft we 



consider the rotating-wave approximation (RWA). 

First, we split the Hamiltonian into time- independent and time-dependent parts: 

A e + Acosu}t TT Tr/ . 
ff = --a x -^ <r z = ff + V(t), (CI) 

Ho = -yffa, - (C2a) 

V(t) = -^f^a, (C2b) 

Here we consider the time-dependent term with a cosine function (in contrast to the sine function used in Eq. @) 
because this choice agrees with the standard convention in the literature. The difference between the two choices is 
simply a shift of wt by 7r/2, and it does not change the time-averaged probabilities which are of most interest to us. 
Then we make a transformation to a rotating frame with the operator W: 



exp f-i J V(t)dt) = exp [i^Y^^j 



W(t) = exp -i / V(t)dt = exp i^a z , (C3) 
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n— — oc 



^4 

r](t) = — smut. (C4) 

This operator matches the wave function in the reference frame ip with the wave function in the rotating frame ip': 
if] = W(t)ip'. Then the Shrodinger equation in the rotating frame reads: iip' = H'ip' with 

H' = W + H W = -^( e ~^a+ + h.c.) - ^a z , (C5) 
a + = | (a x +i(7 y ). Making use of the formula 

oo 

gtesinr _ J n (x) e mT , (C6) 

the new Hamiltonian takes the form 

oo 

H' = - J2 ^f(e-™ t * + + h.c.)-*fa z , (C7) 

TL — — - OO 

A n = AJ n (A/u). (C8) 

Note that the unitary transformation W does not change the level occupation probabilities (the absolute values of 
the spinor components): 

2. Solving the Schrddinger equation in the absence of relaxation 

Let us now study first the system without relaxation. We now look for a solution of the Schrodinger equation 
iip' = H'ip' in the form: 



, /ViexpM^ \ 
V- 1 = , L J\- (CIO) 

1 ^exp(i^) ' 



Actually, this corresponds to the transformation: ip" = W[t)ip' with W(t) = exp(z^jier z ). 

Consider now parameters close to the /c-photon resonance, where ku w Ai5 = |e |. In other words, we now consider 
the frequency detuning 5lj — kui — £o to be small. The RWA consists in the assumption that, in the vicinity of the 
fe-th resonance, fast oscillating terms with n ^ k can be neglected. Then the Schrodinger equation can be readily 
solved. The probability of the upper diabatic state (if the system was initially in the lower diabatic state) is given by: 

-nW/^ l , /.m2 1 A 



p %(t)=\MtW=^ Sua ~+ A 2 (l-cos^/S^ + Altj. (Cll) 

3. Solving the Bloch equations with relaxation 

Now let us take into account rela xation by using the Bloch equations, which include two phenomenological relaxation 
times, Ti and T2 ( Abragaml 1 1 6 lh ■ A convenient parametrization for the density matrix is the following: 



P=U 1 + Z X - iY ). (C12) 

Here we note that the Bloch equations have to be written in the adiabatic (energy) representation, so that the constants 
Ti and T2 describe relaxation to the equilibrium values: Z — > Zq and P21 = X + iY — s- 0. This representation differs 
from the diabatic representation (with the states {(p^, <Pi}) by a small rotation with the angle A/eq <C 1. Here we 
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ignore this difference in the bases and consider the diabatic representation. Then for Z and p2\ — X + iY, the Bloch 
equations with the Hamiltonian H' read: 

Z = -Afelm (e- tk ^p 2 i) - (C13a) 

J-i 

P21 = iA k e ikut Z-ie p 21 -^. (C13b) 

(Again, we have neglected fast oscillating terms with n ^ k.) Here Zq — tanh AM. describes the equilibrium energy 

level distribution at temperature T. After the substitution ^21 = P21 exp(— ikut) = X + iY, we obtain the system of 
equations: 

X = SuY-^, (C14a) 

Y = -SujX + A k Z- — , (C14b) 
T2 

Z = -A k Y-^—^. (C14c) 
J-i 

The stationary solution of these equations is obtained after substituting X = Y = Z = 0. Then the stationary value 
for the probability of the upper diabatic state is given by 

P { £=p 22 = \{l-Z), (C15) 

Z - l + Tl6^ + Tl T2Al Zo - (C16) 
At low temperatures we can substitute Zq m 1 and obtain from Eqs. (IC15IIC16|) the following formula 

Pup = J2 P "P = o zZ ~T , T 2 c 2 j » ( C17 ) 
which is useful for the description of multiphoton resonances at A/e *C 1. 



Appendix D: Floquet theory 



Floquct theory can be applied to any quantum system governed by a time-periodic Hamiltonian. The main idea 
of applying this theoretical approach is as follows: the time-dependent problem is turned into a time-independent 
one at the expense of increasing the effective number of states in the Hilbert space to an infinite number of states. 
Approximations are then made to reduce the problem into a finite-dimensional one, and analytic or numerical results 
are obtained. One advantage of Floquet theory is that it allows for the inclusion of dissipative processes relatively 
easily into the analysis. In this A ppendix we brief ly outline the application of Floquet theory to the problem of LZS 
interf er ometry, closely followin g ( Son et al . 20091). M ore detailed discussions can be found in ( Grifoni and Hanggll . 



2009). 



Il998h . (|Chu and TeFnovl . l2004h . and (|Son et all 

From Floquet's theorem we know that, since the Hamiltonian given by Eqs. ([T]) and 
solutions to the Schrddinger equation with the property that 



is periodic, there exist two 



(Dl) 



and ipj(t) is a periodic function of t with the same period as the Hamiltonian, i.e., 2tt/u. Furthermore, these two 
quasi-periodic solutions are orthogonal and form a complete basis. As a result, knowledge of these states and of the 
quasi-energies e allows one to describe the dynamics for any given initial state. 

Finding the Floquet states and quasi-energies is typically achieved using Fourier analysis. Writing 



(D2) 
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where each £j iTl = (£j, n ,+ > £j.n,-) T is a two-component vector (with components £j, n ,+ and £j.n,-), and writing the 
Hamiltonian (|C1[) as 



A 



.4 



2 2 4 
one obtains the following set of equations governing £j- n : 



a. (e lujt + e~ iut ) 



.4 



-<Tz£ 



A 



j,n-l 



j,n+l- 



(D3) 



(D4) 



The task is now to find £j jra for all values of n (noting that £j >n is composed of the components £j, n ,— and £j,n,+) 
and along with them the quasi-energy 6j (note that the above equations apply to both values of j). If one constructs 

the infinite-dimensional vector £j = {..., £j <n _i t _, ^j in _x,+, £j,n,— > £j,n,+ > £j,n+i,— > £j',n+i,+j •••}> then the above set of 
equations governing £j )n turns into a single eigenvalue problem 



(D5) 



where 



/ • 



-f + (n-l)w 

_ A 
2 
A 
4 







A 

2 

f + (n-l)w 



A 
4 






4 



2 



A 

2 
A 

4 



A 
4 
_ A 
2 

+ no; 


A 
4 






A 
4 





-f + (rc+l)cj 



A 

' 2 







A 
4 
_ A 
2 



f + (n + l)w 



V 



(D6) 



/ 



This eigenvalue problem has an infinite number of solutions, even though we started by looking for only two solutions. 
This issue is resolved by noting that this infinite number of solutions can be divided into two groups, each of which 
contains an infinite number of equivalent solutions; clearly, if there is a solution to the original Floquet problem with 
quasi-energy e, then essentially the same solution can be re- written such that the quasi-energy is given by e + nuj, 
with any integer n. 

The eigenvalue problem derived above cannot be solved in analytic closed form, and approximations have to be 
made. It is common to treat A as a small parameter in a perturbation-theory-like calculation. Taking A = 0, one 
finds the (unperturbed) eigenvectors (|Son et all 120091 ) 



7(0) 





-,n+fc— 1,- 


- Jl 


\2ui 
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-,n+k-l,+ 


= 




& 


— ,n+fc,— 


= Jo 


( — 
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-,n+fc,+ 


= 




& 


-,n+fe+l,- 


= Ji 


V2w 




-,n+fe+l,-(- 


= 





?(0) 



£fc, + .n+fc-l,- 


= 




£k,+,ra+fc— 1,+ 


- 


2ui) 


£fc, + : n+fc,- 


= 




£ft,+,ra+fc,+ 


= Jo(- 


— ) 


£jfe,+,n+k+l,- 


= 






- M- 


±) 

lull 



with quasienergies Ek,± — i £ o/2 — kuj (note that the above eigenvectors are normalized to unity). When two 
quasienergies are almost equal (specifically So,- and £&,+), and assuming that A/lj is small, all others are far off 
resonance from these two, one can make a rotating-wave approximation and construct an effective two-level problem 
that involves the two respective Floquet states. Using the relation 



> J n (x)Jh-n{%) 



n=—co 



(D7) 
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one can calculate the effective coupling strength between the two relevant states. The resulting 2x2 effective 
Hamiltonian is given by 



H = l( so —AJ-k (A/lj) 

F ' RWA 2 ^ -AJ_ fe (A/lj) e - 2ku 



(D8) 



Using this effective Hamiltonian, it is now straightforward to find the Floquet states and quasi-energies. The separation 
between the quasienergies on resonance (i.e., when Sq = kui) is given by \AJk(A/uj)\ and determines the oscillation 
frequency between the two states of the TLS, in clear analogy to the RWA calculation of Appendix C. One can 
proceed wi t h the calculation in order to obtain more accurate approximations to the Floquet states and quasienergies 
( Son et al ., 2009). but we shall not do this here. 



In this appendix we have re-derived results that we had already derived using the adiabatic-impulse model and 
the RWA. The use of Floquet theory in this context can therefore be seen as just an alternative approach that can 
be used to obtain results that can be obtained using other methods. One advantage that the Floquet approach 
has over the adiabatic-impulse model and the RW A is that decoherence can be introduced into the problem in the 
Floquet approach following a standard formalism ( Grifoni and Hanggl Il998l ). whereas the Bloch equations that we 



have used in Sec. II can be seen as a phenomenological approach to including decoherence in the theoretical analysis. 
The question of the effects of decoherence on LZS interferometry in the Floquet approach is currently under study 
(|Hausinger and Grifoni l201flh . 



Appendix E: Dressed-state picture: quantized driving field 

In all our discussion above, we have treated the driving field as an externally applied, classical driving signal. An 
alternative theoretical approach that can be used to study the problem of a driven quantum system treats the driving 
field quantum mechanically. The quantum system under study is then enlarged from simply a two-level system to a 
two- level system and a harmonic oscillator that can have any number of excitations (i.e., photons). As in the Floquet- 
theory approach, a major advantage that justifies the expansion of the Hilbert space is the simplification associated 
with turning the time-dependent Hamiltonian into a time-independent one. In fact, the algebra encountered in the 
Floquet-theory calculation is identical to the one encountered in the dressed-state analysis. The difference between 
these two approaches is therefore simply the physical interpretation of the mechanisms at play during the system 
dynamics. 

The equivalence between the Floquet approach and the dressed-state approach can be seen by considering the 
Hamiltonian of the dressed-state picture: 



H 



— a z + wo 1 a - 



.4 



(EI) 



where a and a) are the annihilation and creation operators for photons in the electromagnetic field. When written in 
matrix form, the above Hamiltonian is giv en by Eq. ( D6I) d e rived in the Floqu et approach. More detailed studies of 
the dressed-state picture can be found in [ dLiu et aLl . l2006h . dGreenberd. l2007ri . dWilson et aZ.1 . 120071 12010T) ]. 
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FIG. 11 (Color online) Experimentally realized L andau- Zener-Stiick elberg (LZS) interferometry. The p anels from top to bo ttom 
present the resu l ts of the following articles: (a) (|Qliver et aU 120051 ), (b) (|Sillanpaa et all 120061 ), (c) l| Wilson et all 120071 ), (d) 
(jlzmalkov et aLll2008h . Schematic diagrams of the circuits used are shown to the left, while results for the LZS interferometry 
are presented to the right. A more detailed description of th e experiment s can be found in the main text and, of course, in 
the respective original articles. Figur e (a) is reprinted from IjOliver et al. 1 120051 1 with permission from AAAS. Figure (b) is 
repri nted from ( Sillanpa a, et all 120061 ) with permission; copyright (2006) by APS. Figure (c) is rep rinted from (|Wilson et all 
2007) with permission; copyright (2007) by APS. Figure (d) is reprinted from (llzmalkov et all [20081 ) with permission; copyright 
(2008) by APS. 



